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Abstract: 

We propose a simple parametrization for the form factors in heavy to light decays which 
satisfies the heavy quark scaling laws and avoids introducing the explicit "dipole" form. 
The illustration is provided on the set of lattice data for B — > 7r semileptonic decay. The 
resulting shape of the form factor also agrees with the QCD sum rule predictions. 
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1. Preliminaries A lack of the precise information about the shapes of various form 
factors from the first principles of QCD was, and still is the main source of uncertainties 
in extraction of the CKM parameters from the experimentally measured exclusive decay 
modes |] . This problem is particularly pronounced in the heavy — > light processes because 
the kinematically accessible region is very large. In such a situation it is not only important 
to know the absolute value of a particular form factor at one specific point (traditionally 
at q 2 = 0), but also its behaviour in q 2 , which (when suitably integrated over the whole 
phase space) gives the physically measurable branching ratio. So far, there were very many 
studies devoted to that issue. In spite of the recent progress, there is still no method based 
on QCD only, which could be used to describe the complicated nonperturbative dynamics in 
the whole physical regionQ Most of the approaches agree that the functional dependence of 
the heavy — > light form factors is 'distorted' in the region of not so high q 2 (or equivalently, 
for moderately and very large transfers, q 2 ). In other words, it is different from what one 
would obtain by invoking the nearest pole dominance. Although somewhat intriguing, this 
discrepancy with the pole dominance hypothesis is however expected: the low q 2 region is 
very far from the first pole, so that many singularities (higher excitations and multiparticle 
states) become 'visible' too, i.e. their contribution becomes sensible (comparable to that 
of the first pole). One would obviously like to quantify these effects of higher states. 

For definiteness, we consider the form factors f +,0 {q 2 ) which parametrize B — > iriv 
decay as 



(n(p')\V»(0)\B(p)) = (p + p'-q ml q2 ml )^f + (q 2 ) + ^-^q»f°(q 2 



where = wj^b. f + (q 2 ) is the dominant form factor, i.e. the only one needed for the 
decay rate 



/7F IT/ I 2 



in the case of a massless leptonQ (& = e, /i). The subdominant form factor f°(q 2 ), measures 
the divergence of the vector current and its contribution is proportional to m 2 in the decay 
rate, thus it is irrelevant for present experiments Its knowledge is however important 
in lattice analyses, as well as in the approaches relying on the soft pion theorem (in which 
the form factors are calculated at large q 2 ), because it can help to constrain the dominant 
f + (q 2 ), through the kinematic condition 

/ + (0) = /» . (3) 



2. Constraints Beside this important constraint (0), in the heavy quark limit {uib — * oo) 
and near the zero-recoil point {q 2 ~ 0, or equivalently q 2 ~ q^ax) where the pion is soft, 

1 The recent reviews about the form factors in heavy to light decays can be found in Refs. |^|, |j. They 
also contain an exhaustive list of references. 

2 X(t) = (t + m% — to 2 .) 2 — 4m 2 3 m 2 , is the usual triangle function. 
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one encounters the well known Isgur-Wise scaling law j|]f] 



(7r(p')\V^\B(p)) oc ^ ^ 1 + 0(l/ro fl ) 
which in terms of form factors reads: 

f + (<l 2 - ?max, m B) ~ > A*? 2 ^ <?Lx> ™b) ~ ^= ■ ( 4 ) 

As first noticed in Ref. ||, one may use the low energy theorem to get one more constraint, 
often invoked in Kg& decays. Namely, in the soft pion limit (p' — > and m 2 — > 0), the 
r.h.s. of ([ip simply becomes f°{p 2 )p IJ ,, while the l.h.s. from the LSZ reduction of the pion 
amounts to (fn/ fn) Pp, so that one deduces: 

/VI) = t-- (5) 

Jit 

Next scaling law comes from the asymptotic limit of the light cone QCD sum rules which 
are applicable in the low q 2 region. It was first shown in Ref. 0, that the form factors in 
heavy to light decays satisfy 

f +0 (q 2 ^0,m B )~m B 3/2 . (6) 

This scaling law was recently derived also in the framework of the large energy effective 
theory 0. @ will provide us an important additional constraint in what follows. 

Our knowledge of the form factors goes a little bit beyond. We know that the Lorentz 
invariant form factors are analytic functions satisfying the dispersion relations: 

W) = ir dt w^L , (7) 

it Jto t — q A — le 

Res g2=w 2 fl J+(g 2 ) i ,oo Im/+(t) 



f+tf) = — q : + - / dt — v^ 2 - • (8) 

J V ; m 2 B * - q 2 tt Jto t - q 2 - te w 

The imaginary parts are consisting of all single and multiparticle states with the same 
quantum numbers as f + (t) (f°(t)), i.e. J p = (0 + ), thus a multitude of poles and cuts 
situated above the Bit production threshold, t = {tub + m-ir) 2 - However, below t (but 
above the semileptonic B — > tt region) there is a B* vector meson (rag. = 5.325 GeV) 
contributing to the f + {t) form factor. The residue at that pole is the product of the B*Btt 
coupling and the coupling of B* to the vector current: 



1 

l b* j yi ' 2 



Res g2=m 2 D J + {q 2 ) = -m B *f B *gB*Bn , (9) 



3 The overall logarithmic correction proportional to (a s (rrih) / a s (fi)) (arising from the perturbative 
matching between the QCD and the heavy quark vector current) are irrelevant for our discussion. 
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where the standard definitions were used: 



(0|VgB*(p,e)> = /fl.mB.e, , 
(B-{p)ir + {q)\B* Q {p + q,e)} =g B * B «{q-e) . (10) 

While the lattice is providing us with better and better estimate for /#* ||, the value of 
9b* Btt remains vague. From the above definitions, one immediately finds that Res g 2 =m | t f + (q 

scales as m 3 B 2 . 

The simple pole dominance would mean that only the first term in the r.h.s. of (^) 
survives, whereas all the other states which can couple to the vector current and are 
above to, would eventually cancel. Moreover, one would encounter the discrepancy when 
trying to reconcile the pole dominance ansatz with (|6|), because the pole dominance gives 
/+(0)~m B 1/2 . 

The easiest way out was to suppose the pole behavior for f°(q 2 ), and the double pole 
(dipole) one for f + (q 2 ) |, [K|: 

/V) = i , /V) = n /( 2 °) i nr - (ID 

1 — q /m Q+ ( 1 - q l jm\_ ) 

This choice basically satisfies all constraints but it is not informative at all: we know that 
we have a pole at B* (and more importantly, we know its position) which we trade for 
a better fit with data, without getting any supplementary physical information. For that 
reason, we wanted to propose a new parametrization which would satisfy all the constraints 
mentioned above, and use the fact that we know where the first pole is. In this way the fit 
with data will turn to be far more informative. 
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3. Proposal One should always start from the fact that f + {q 2 ) does have a pole at m B *. 
The contribution of other higher states (for which we do not know positions, nor residues), 
can be parametrized as an effective pole whose position is to be determined relatively to 
the first pole. In other words, we propose: 

where we introduced x = q 2 /m B *, and parameters a and 7 - positive constants which scale 
with tub- Obviously, cbvp? b * is the residue of the form factor at B*. Note that / ym 2 B , is 
the squared mass of an effective 1~ excited _B*'-state, which is equal to (mg. + A) 2 ~ 
m 2 B * + 2Am#», where A does not depend on the heavy quark in the heavy quark limit. 
Therefore, 7 ~ 1 + 2A/ms». We now check the scaling laws. At q^axi one nas: 

(o?7 \ 
1 + Ti \ (I 3 ) 
(l-y)m B J 
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Since c B scales as ~ m^ 1 ' 2 , then in order to fulfill the heavy quark scaling law (f|) 
f + (%max) ~ m B 2 ) 1 — 7 can at most scale as l/m_B. That is consistent with the above 
observation on 7. It is also clear that in this counting, aj scales as a constant. On the 
other hand, at q 2 = we have: 

/ + (0)=c B (l-a)~m B 1/2 (l-a)~ m 5 3/2 , (14) 

which means that (1 — a) scales as ~ 1/itlb- This shows why we have chosen the sign 
"— " between the two terms in (|12|) . In fact, it is easy to see that both a and 7 scale as 
1 + const. /m B + ■ • •. To build /°, we take into account /°(0) = / + (0). For convenience, 
we can measure the position of the effective pole to f° in terms of x = q 2 /m B *, and write 



/•(,) = -22(!z°L . (15) 

1 - T3 



One can easily see that f°(qmax) ~ m B^ as ^ should, while /°(0) ~ °m B ^ 2 by construction. 
Using the same argument as for 7, it is obvious that (3 also scales like 1 + 0(l/m B ). 

These four parameters (c B , a, (3, 7) contain a considerable information on the process: 
cb gives a residue of the form factor at the pole B*, a measures the contribution of the 
higher states which are parameterized by an effective pole at m B * = j-m B * . Contributions 

cff 

to the scalar form factor are parameterized by another effective pole whose position is at 
m%>» = ■ m%*. The situation with the form factor f + resembles the situation with the 

B cff r ° J 

electromagnetic form factor of a nucleon, where the simplest p-dominance is not sufficient 
to describe a dipole behavior of the form factor at large q 2 , and one needs to introduce at 
least one extra excited vector meson. 

Present lattice analyses have limited statistical accuracy and accessible kinematical re- 
gion for the form factors and for a description of lattice data a more constrained parametriza- 
tion is required. In the limit of large energy of the light meson (large recoils), in the rest 
frame of the heavy meson 

771 n / q 2 m 2 \ m B ^oo Uln , , . 

E=^-[l-\ + ^) ~ 1-x . 16 

2 \ m B m B J 2 

In Ref . , it was shown that in the heavy quark limit and for large recoils (m B — ► 00 and 
E 00), there exists one nontrivial relation between the two form factors, namely: 

f° = —f + . (IT) 
m B 

In this limit, and by using the above counting argument with a = 1 — ao/rriB, 7 = 
1 + •yo/rriB, and (3 = 1 + Po/n^B, we have 



+ _ c a 
" 2E^f^ 
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(18) 



Therefore, to satisfy the relation fll7D it is necessary to have 70 = ao- It is reasonable then 
to take, a + 7 = 2, but also a = I/7, which satisfy ( fPTD up to terms 0{l/m 2 B ). In what 
follows, we use = 1/7. 

To summarize, the recipe is the following one: in lattice analyses, one first calculates 
rriB* (from the corresponding two-point correlation function)^ to measure x, and then fits 
the data according to: 



1 




1 


(-})■ 


f°(x) 


CB 


;i-a) 


1 




1 


(1 — x) (1 — ax) , 


f + (x) 


cb 


I -a) 



(19) 



which makes altogether three parameters (eg, a, /3) to perform this fit. 



3. Illustration Now, we would like to show how it works on the specific example with 
the set of lattice data. Before the new (systematically improved) results are released^, we 
decided to use the most recent available results, those of JLQCD JTT|. Since our purpose 



q 2 [GeV 2 ] 


20.7 


19.2 


18.3 


17.4 


16.4 


15.4 


X 


0.70 


0.65 


0.62 


0.59 


0.55 


0.52 




0.72(3) 


0.71(3) 


0.67(2) 


0.66(4) 


0.65(4) 


0.60(7) 




2.09(8) 


1.52(8) 


1.41(7) 


1.28(10) 


1.20(12) 


1.11(17) 



Table 1: Lattice results obtained by the Hiroshima group, JLQCD ||1 1|| . 



is only to illustrate the parametrization (|i9|) , we will not further comment on the technical 
details used in the calculation of JLQCD. However, one should note that the lattice groups 
reported results by assuming the negligible SU(3) breaking effects, which is also supported 



by the QCD sum rule calculation [| [12"], as well as by experiment with D-mesons [JTJJ]. That 

In experiment, one fixes it to the physical mass of the vector meson B*. 
5 APE and UKQCD groups, are already analyzing their new respective data. 

6 Note that the LCQCD sum rules predict larger SU(3) breaking which is the consequence of the 
normalization of the light meson wave functions (implying that fg^ K (0)//s_ ff (0) ~ fx/fir)- This is in 
conflict with the findings of lattice calculations 0. Further research clarifying that issue is needed. 
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means that the pion used in the JLQCD calculation is consisted of degenerate quarks of 
mass near the strange one. Also, the spectator quark in B meson is the V quarkQ. 

First, we fix the mass of mg. = 5.46 GeV, as obtained in the same calculation of 
JLQCD. In realistic situation, one of course takes the experimental mg. = 5.32 GeV. 
Then we fit the data with the expressions ([19]), i.e. by minimizing: 

AVn - (i - f ) v - *)]), - (i - «* 



x 2 =n + 



ft: 



A S(l/f% ) V <W[/+(1 -*)]), 

from which we obtain the following result: 

/+'°(0) = 0.38(8) , ft = 0.54(17) , j3 = 1.4(9) . (20) 



The best fit and the lattice data are depicted in Fig.|I|. These conservative errors in ( ^0|) 
are estimated from the correlated fit, corresponding to A% 2 = 3.53 (68% CL.) W%, which 



defines the ellipsoid of errors (three parameters). The errors we quote in (|20|) correspond 
to the extrema of that ellipsoid. As we already mentioned, by our parametrization, besides 
the first pole in f + (x) all the other singularities are replaced by an effective pole, which 
is situated at mg. = m B *^j = 1.3 jlal tub*- The effective pole for f°(x) is aX m B * = 
m B *VP = 1.2±° fm B *. 

An important advantage of this parametrization is its transparent physical meaning. 
For example, after identifying /(0) = c B (l — ft), we may also extract the information on 
the residue of the form factor / + at B* , i.e. about the coupling g B * Bn : 

c B =0.8(3) => g B , Bw = 1.6(6)™*.//*' • (21) 

If we take the physical value m B * = 5.32 GeV, and f B * ~ 0.2 GeV, we see that the resulting 
9b* Bit — 42 ± 16, which is compatible with most of today's estimates of this coupling (see 
Tab. 2 in Ref. A statistical (and more importantly, a systematic) improvement of the 
lattice data will provide more reliable estimate of this coupling. 

The ellipsoid of errors is of the cigar shape so that one can reduce the errors and 
consequently get some additional information, by making reasonable assumptions. One 
can for instance make the following exercise. We saw that the position of the effective pole 
of the f + (q 2 ) function is at 1.3±g;| m B *. If we now fix it to (1.15, 1.20, 1.25, 1.30, 1.35) xm B *, 
we obtain the results listed in Tab. 2. We notice that the data tend towards larger values of 
the mass of the scalar effective pole as compared to the expected mass of the lowest scalar 
meson Q, /3 — 1.17, though they are consistent with this value. Note that an expected mass 
of the first excited vector state according to theoretical estimates should be close to 6 GeV, 
which again corresponds to a smaller value of 7 ~ 1.3 (or a « 0.7), but within the errors 
agrees with the values found in our analysis. 

7 To our knowledge, among recent studies, only the APE collaboration made the extrapolation to 
the chiral limit. 

8 The scalar (0 + ) meson is expected in the region of the measured bump [|l6| of B** states: m^.. = 
5.73 GeV. 
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m B * Jm B * 

en ' 


1.15 


1.20 


1.25 


1.30 


1.35 


a 


0.76 


0.69 


0.64 


0.59 


0.55 




n r KC\(~\ \ 

U.OU^l J 






U.OU^l J 


n ^Rfi ^ 

u.oo^ i j 


9B* Bit /b* / m B* 


2.45(5) 


2.10(4) 


1.90(4) 


1.76(4) 


1.65(4) 


a 


1.12(6) 


0.19(7) 


1.26(8) 


1.32(9) 


1.39(10) 


m B * Jm B . 


1.06(3) 


1.09(3) 


1.12(4) 


1.15(4) 


1.18(4) 


}b/ fn 


2.5 ± 1.0 


1.8 ± 0.6 


1.6 ± 0.4 


1.4 ± 0.3 


1.3 ± 0.3 



Table 2: Results of fits for a fixed 7. The last row is obtained using Eqs. @ and fllSD in 
which we took x = 0.983. 

Now, when we have fixed the parameters by the lattice data (pTj), we would like to 
compare the resulting shape of the form factor f + (q 2 ) to the result of the light cone QCD 
sum rule, which is applicable in the low q 2 region. To do that we define f + (q 2 )/ f + (0), 
and plot them both in Fig. |2|. We see that our parametrization, constrained by the lattice 
data at large g 2 's, indeed reproduces the shape of the form factor at low g 2 's as predicted 
by QCD sum rule. To be consistent, in the comparison we used the physical m B * — 5.32 
GeV. 

4. Summary In this letter we proposed the physically more informative way to fit the 
data for the form factors in heavy to light decays. This is illustrated on the example 
of B — ► 7r semileptonic decay, and can be also used in D — > tt, as well as in the other 
heavy — ► light transitions. Initially, the parametrization for two form factors contains four 
parameters (|P2"D , fll^D , which we further reduced to three ((19]) . All parameters scale with 
the heavy mass as const. + 0(l/m B ), except for c B . The proposed parametrization gives 
us the information on the singularities which decisively influence the shape of the form 
factor in the semileptonic region, and allows us to extract the value of the residue at the 
nearest pole (i.e. to determine the value of ^-coupling). 

On the specific example, we have shown that when fixing the parameters in our parametriza- 
tion by the lattice data in the large q 2 region, one reproduces (to a good precision) the 
shape of the form factor obtained by QCD sum rules at low and intermediate g 2 's. 

We believe that this parametrization may be useful in the forthcoming analyses of 
systematically improved lattice data, as well as in experimental studies of exclusive heavy 
— > light decays. 
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Figure 1: Fitting the lattice data using the parametrization fll^p. Note that the fit of f + 
form factor is constrained by the precise data for f°. For easiness, only the central curves 
(without errors in parameters) are displayed. 




Figure 2: The shape of the form factor from our parametrization ([191), constrained by the 
lattice data at large q 2 , is compared to the QCD sum rule predictions at "low" q 2 ( Eq. (132) 
in Ref. || for the dashed curve, and Eq. (23) in Ref. for the dotted curve) . A very good 
agreement is found. 
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